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Abstract: We monitored the Ca II H and K lines of HD 179949, a notable star in the southern 
hemisphere, to observe and confirm previously identified planet induced emission (PIE) as an effect of 
star-planet interaction. We obtained high resolution spectra {R ~ 53, 000) with a signal-to-noise ratio 
S/N > 50 in the Ca II H and K cores during 10 nights of observation at the McDonald Observatory. 
Wide band echelle spectra were taken using the 2.7 m telescope. Detailed statistical analysis of Ca II 
K revealed fluctuations in the Ca II K core attr ibutable to planet in duced chromospheric emission. 
This result is consistent with previous studies by Sh kolnik et al.l l|2003l ). Additionally, we were able to 
confirm the reality and temporal evolution of the phase shift of the maximum of star-planet interaction 
previously found. However, no identifiable fluctuations were detected in the Ca II H core. The Al I 
A3944 A line was also monitored to gauge if the expected activity enhancements are confined to the 
chromospheric layer. Our observations revealed some variability, which is apparently unassociated with 
planet induced activity. 
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1 Introduction 

Observational results indicate that the large major- 
ity of planets outside t he solar system is hosted by 
main-sequence stars (e.g. lButler et al1l2006f l. which are 
known to possess an extended outer atmosphere con- 
sisting of a chromosphere, transition region and corona 
as well as a dynam ic stellar wind region (e.g. iLinskvl 
Il980l : ISimo"iilll986l) associated with comp lex magnetic 
structure (e.s. ISchriiver fc Zwaan|[2000l ). In various 
systems, consisting of a star and a close-in extrasolar 
giant planet (CEGP), interaction between the planet 
and stellar atmospheric structure has been identified 
that is often broadly classified as "planet-induced (stel- 
lar) emission" (PIE). The reality of thi s phenomenon 
has been proposed bv lCuntz et all (|200G| ') who also pro- 
vided a detailed ranking of the relative strength of this 
effect for the var ious star-planet sys tems known at the 
time. Moreover, ICuntz et all (l200Gf ) explored the pos- 
sibility of tidal interaction and magnetic interaction 
with the latter distinguished from the former by two 
activity maxima per stellar circumference instead of 
one. 

The PIE effect was subsequently discovered bv lShkolnik 



(12003) in regard to the HD 179949 system using high- 
resolution {R = 110, 000) and high signal-to-noise ratio 
spectra based on three observing runs at the Canada- 
France-Hawaii Telescope (CFHT). They also concluded 
that the PIE effect is expected to be magnetic in nature 
rather than caused by gravitational sta r-planet inter- 
action . Subseque nt work was pursued by lShkolnik et all 
(120051 ) ■ Shkolnik eTal] (|2008t l and others, resulting in 
the identification of the PIE effect in at least six sys- 
tems, which are: HD 209458, v And, r Boo, HD 179949, 



HD 189733, and HD 73526. 

The latter study also provides evidence of the on/off 
nature of the PIE effect, also previo usly found in the 
V And system l|Shkolnik et al.l [20081 ). which according 
to the authors is likely attributable to the changing 
stellar magnetic field structure thr oughout the stel- 
lar activity cycle (e.g. iLanzal [2OIOI ). An alternative 
or supplementary explanation of the significant time- 
dependent variations of the PIE phenome na refers to 
the flare- type nature of th e interaction (e.g. lSaar et al.l 
l2004l : iMcIvor et alj|200^ . Recently, a detailed time- 
dependent MHD simulation for the interaction of the 
stellar magnetic field and wind with the planetary mag- 
netosphere and outflow was given bv lCohen et"an (|201ll ) 
for HD 189733 as an example. Their results show re- 
connection events occurring at speciflc planetary or- 
bital phases, causing mass loss from the planetary mag- 
netosphere able to generate hot spots on the stellar sur- 
face. The simulations also demonstrate that the sys- 
tem has sufficient energy for the hot spots to be visible 
in Ca II lines, which are expected to show stochastic 
behavior. 

There were also various attempts to identify the 

^Hli .[ effect in other spectral regi mes, including t he ra - 
dio regime, infrared and X-rays. iBastian et all (|2000l 'l 



conducted a search for radio emission for six extra- 
solar star-planet or brown dwarf systems that included 
four systems with the planet or brown dwarf closer 
than 0.1 AU from its host star. No detections were 
made, however, which was attributed to physical or in- 
strumental reasons. Subsequent work mostly focused 
on the derivation of upper limits of eme rgent radio 
emissio n from known syste ms given by. e.g. [Farrell et al] 
(|200^ '). IZarkal (120071 ) . and iGriefimeier et all ^20ff^ . 
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Tabic 1: Stellar and Planetary Parameters 



Parameter 


tiU 179949 


Keierencc 


Spectral Type 


F8.5 V 


a 


RA 


19*^ 15™ 33.2278'' 


b 


DEC 


- 24° 10' 45.668" 


b 


V 


6.25 ± ... 


b 


Tcff (K) 


6202 ± 52 


Ribas at al. f2003) 


(Mq) 


1.24 ± 0.10 


Tinnev et al. l'2001~) 


Distance (pc) 


27.0 ± 0.5 


b 


Age (Gyr) 


2.05 ± ... 


Donahue C 1993) 


R* (Rq) 


1.193 ± 0.030 


Ribas et al. f2003~) 


[Fe/H] 


0.226 ± 0.050 


Gonzalez & Laws ( 2007) 


ap (AU) 


0.0443 ± 0.0026 


Butler et al. (2006) 




0.022 ± 0.015 


Butler et al. (2006) 


Porb (days) 


3.092514 ± 0.000032 


Butler et al. (2006) 


MpS'mi (Mj) 


0.916 ± 0.076 


Butler et al. (2006) 



"data from SIMBAD; see |h ttp : / /simbad . u-strasbg . f rl 

adopted from the Hipparcos catalogue; see lESAI (|1997l) . 



Anot her effort to identify t he PIE effect was under- 
taken bv lSaar &: Cunt3 (|200ll 'l who executed a search 
for periodicities in the Ca II infrared triplet emission 
for a sample of stars with close-in giant planets; how- 
ever, no unequivocal identification was made. Sub- 
sequently, an effort to search for the PIE effect in 
the UV using the Goddard High Resolution Spectro- 
graph (GHRS) on board of the Hubble Space Tele- 
scope (HST) was rendered impossible as the GHRS 
was re moved during H ST Servicing Mission 2. Re- 
cently, iGonzale^ (|201ll ) reported differences regarding 
the chromospheric activity between samples of stars 
with and without planets, although the interpretation 
of these trends, if confirmed, remains unclear. 

iKashvap et al.l (l2008l ) presented a statistical sur- 
vey of the X-ray fiuxes from stars with close-in planets; 
arguably, this survey has been highly useful as these 
fiuxes were found enhanced by 30% to 40% on average 
over typical fiuxes from similar stars with planets that 
are not close-in. On the other hand, contrary evidence 
was subsequently identified bv | Poppenhaeger et al.l (120101) , 
iPoppenhaeger fc Schmit3l|201ll ). and lCanto Martins et al.l 
(|201ll ) ; nonetheless, these studies, if found conclusive, 
do not exclude the possibility of enhanced X-ray fluxes 
for stars with planets that are very close-in (i.e., < 
0.05 AU). Specifically, d irect X-ray observations of the 
HD 179949 system by ISaar et al.1 (,2008) revealed a 
~30% increase in the X-ray flux over quiescent levels 
coincident with the phase of the Ca II enhancements. 
Additionally, the observations indicate a trend for the 
emission to be hotter at increased fluxes, conflrmed by 
modeling, showing the enhancement at ~1 keV com- 
pared to ~0.4 keV for the background star. 

In our present study, we will again focus on the 
HD 179949 system. HD 179949 is an F8.5 V star 
at a distance of 27.0 ± 0.5 pc located in the south- 
ern hemisphere in the constellation Sagittarius with 
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Figure 1: Orbital position of the planet around 
its host star (HD 179949) during the observational 
run. Note that 4> — corresponds to the planetary 
position between the host star and Earth (sub- 
planetary point). The numbers indicate the dates 
in April 2006 (local time) when the observations 
were obtained. 
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Figure 2: Single flat-fielded and blazed spectrum 
of HD 179949 taken on April 5, 2006 with an ex- 
posure time of 1200 s. We depict the Ca II regime 
(3915 ~ 3985 A) of the wide band echelle spectrum 
{top) and the associated S/N ratios {bottom). 

a magnitu de of V = 6.25. I ts mass is estimated as 
1.24 Mq (jTinnev et all bOOlh based on the interpo- 
lation of evolutionary t racks previously computed by 
iFuhrmann et all (|l99lll998h . The age of HD 179949 
has been estimated as 2.05 Gyr (see Table 1) based 
on the strength of the c hromospheric component of 
the C a II emission, -Rjjk l|Noves et al.l [l984: Donahue 
119931 '). Note that HD 179949 is a solar-type star of 
moderate age, implying that its inherent chromospheric 
emission and variability has considerably subsidecQ, 
thus adding to the credibility of any PIE determina- 
tions. By employin g the Doppler velocity technique, 
iTinnev et al.l (120011 ) discovered as part of the Anglo- 
Australian Planet Search a "Hot Jupiter" hosted by 
HD 179949 located at a distance of about 0.044 AU 
(7.99 stellar radii). The planetary orbit has a very 
low eccentricity that is essentially consistent with the 
circular case (see Table 1), allowing us to ignore the or- 
bital ellipticity, if any, in the following. The planetary 
mass parameter is given as Mp sini = 0.916 ± 0.076 
Mj. 

Significant evidence that the activity enhancement 

^See Guinan ct al. (2003,) and iGiidell ll2007l 'l. among 
other contributions to the "Sun in Time" project, for stud- 
ies on the decline of chromospheric emission on stellar evo- 
lutionary time scales. 



in HD 179949 is planet-induced is based on the finding 
that the stellar rotation period Prot almost certainly 
disagrees with the orbital period of the planet. There- 
fore, it would be difficult to explain emission on the 
star that essentially "follows" the CEGP without a 
planet-related cause. Unfortunately, HD 179949 does 
not have a robustly measured Prot fsee'Wolf fc H armanec I 
I2OO4). However, following Saar et al. (2004) and refer- 
ences therein, it can be argued based on the strength of 
the Ca II HK emission, the mean Mount Wilson index 
< S > as well as stellar evolution constraints that Prot 
lies between 7 and 10 d, thus being distinctly different 
from Porb ^ 3.092 d (see Table 1). 

The PIE effect has various interesting astrophysi- 
cal applications as it, in principle, may allow the iden- 
tification and characterization of magnetic fields and 
dynamo activity in exo solar Jupite r-type planets (e.g. 
Id ct al. 2004; Gricfimci eTet al.l200 4: Erkaev et al. 20051: 
Zarka 2007; Lanza 2008; Coh en et al.ll2009t : lScharfl2o{or 
Lanza 2011; Cohen ct al. 20^ In addition, it may 
also be useful for probing the inner region of th e stel- 
lar winds (e.g. lPreusse et al.|[2005l : lMcIvor et al.|[2006. ') . 
This is a further reason to continue monitoring the PIE 
effects in different spectral regimes and by using differ- 
ent type of instruments including the 2.7 m Harlan J. 
Smith Telescope at the McDonald Observatory. The 
latter allowed us to continue studying the Ca II H and 
K line cores in HD 179949 and to compare our findings 
with the previous results by Shkolink et al. (2003) and 
subsequent work. Our paper is structured as follows: 
In §2, we describe the employed methods including the 
data reduction. Our results and discussion are given 
in §3. In §4, we present our conclusions. 

2 Methods and Data Reduc- 
tion 

In order to study the PIE effect in the HD 179949 sys- 
tem, we obtained high resolution {R = A/AA ~53,000) 
spectra using the McDonald Observatory's echelle spec- 
trograph mounted on the TK3 light path of the Coude 
focus of the 2.7 m (107 inch) Harlan J. Smith Tele- 
scope. We obtained observations for six (almost) con- 
secutive nights in April 2006, which allowed us to ac- 
quire data points spreading over nearly three planetary 
orbits (see Fig. 1 and Table 2 for details). Full optical 
bandwidth (3400 - 10,900 A) were recorded in at least 
67 orders by a Tektronix (2048 x 2048), CCD cam- 
era. The Ca II H and K lines were carefully centered 
in order 56. Furthermore, Th-Ar spectra were taken 
frequently (i.e., typically in less than 2 h intervals) to 
allow accurate wavelength calibrations. The signal-to- 
noise (S/N) ratio of the spectra exceeded 200 for the 
continuum, and it was close to 50 in the Ca II K core 
(3933 A), and about 80 in the Ca II H core (3968 A) 
(see Fig. 2). 

The data reduction was performed by using the 
IRAF standard packages IMRED, CCDRED and ECHELLE. 
The data reduction procedures can be summarized as 
follows: Combined biases were subtracted from all stel- 
lar, arc, and fiat-field exposures to remove the baseline 
noise. Scattered light was removed from all stellar and 
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Table 2: Details of Observations 



Number 


Date 


Time 


Exposure Time 


S/N 


Seeing 






(UT) 


(s) 




(arcs) 


1 


April 5, 2006 


11:31 


1200 


53 


2.6 


2 


April 5, 2006 


11:54 


1200 


56 


2.6 


3 


April 6, 2006 


11:17 


1200 


43 


4.1 


4 


April 6, 2006 


11:40 


1200 


47 


4.1 


5 


April 7, 2006 


11:30 


1200 


37 


5.0 


6 


April 7, 2006 


11:53 


1200 


45 


5.0 


7 


April 8, 2006 


11:13 


1200 


50 


2.3 


8 


April 8, 2006 


11:36 


1200 


51 


2.3 


9 


April 8, 2006 


11:59 


1200 


70 


2.3 


10 


April 9, 2006 


11:59 


1200 


44 


2.1 


11 


April 11, 2006 


11:31 


1200 


51 


1.7 


12 


April 11, 2006 


11:54 


1200 


66 


1.7 



Table 3: Ephemerides of HD 179949b 



Date 


JD 


Orbital Period 


Reference 






(d) 




07 / 1998 


2451001.510 


3.092514 


Butler et al. f2006)° 


09 / 2000 


2451793.980 


3.0925 


b 


07 / 2002 


2452479.823 


3.09285 


Shkolnik et al. C2003) 


09 / 2003 


2452894.110 


3.09246 


Shkolnik et al. C2003) 



"based on data obtained bv lTinnev et al. 1 (l200lh . 
'^^http : //www ■ exoplanet . eu 
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Figure 3: Residual flux of the Ca II K core of 
HD 179949 between April 5, 2006 and April 11, 
2006. The residual flux is computed from the av- 
erage mean and smoothed over 21 pixels. The re- 
sults are indicated by solid (April 5), dotted (April 
6), dashed (April 7), dash-dotted (April 8), and 
long-dashed (April 11) lines. Note that the data 
of April 9 have been disregarded. 



flat field exposures. All stellar and arc exposures were 
divided by combined flat fleld for normalization of the 
sensitivity levels of the CCD chip. Spectra were ex- 
tracted from the 2D images to 2D spectra. Heliocen- 
tric correction was applied to each spectrum to express 
the observation time considering that the observer's 
location does not coincide with the center of the Sun. 
Differential radial velocity corrections due to the rota- 
tion and orbital velocity of Earth were applied to each 
spectrum. The combined flat field was processed like 
a stellar image to obtain blaze functions for each ob- 
servation night. Blaze functions were used to remove 
from each final spectrum the low-order curvature as a 
result of the intensity variation along the orders. 

Each reduced spectrum is carefully checked for po- 
tential contaminations. Trace amounts of sunlight con- 
tamination are detected in the last exposures of each 
observing night as HD 179949 is observed at low alti- 
tudes on the Eastern horizon just a few hours before 
dawn. These exposures were discarded to avoid any so- 
lar contamination of our spectra. Although the night 
of April 9 originally indicated a critical phase to de- 
tect star-planet interaction, this entire data set had to 
be discarded as the cloud coverage caused large noise 
level in the spectra. 

Ne xt we followed the procedure of IShkolnik et all 
(|2003l ) previously used to unveil the PIE effect for the 
HD 179949 system from optical spectra. The follow- 
ing steps were undertaken: (1) From each spectrum 
increments of 7 A centered at the Ca II H (3968 A) 
and K (3933 A) cores, respectively, were extracted. 
(2) To normalize each spectrum, the end points of 
each extracted spect ra were set to 1. As stated by 
IShkolnik et al.l (|2003l ). this is the best way to normal- 
ize the spectra in the wavelength regimes, such as H 
and K lines, where there is no clear continuum. (3) 
The normalized spectra were grouped by date and co- 



Figure 4: The solid line indicates the Mean Ab- 
solute Deviation (MAD) of the Ca II K core of 
HD 179949 as observe d. The MAD va l ues ar e cal- 
culated as described bv lShkolnik et al.l ( 20031) . not- 
ing that MAD — IZiLi \datai — mean\ based 
on N spectra. The unit is the relative intensity 
attained as a fraction of the normalized flux. The 
dashed line represents the mean spectrum of the 
Ca II K core (scaled accordingly) to show that the 
activity is confined, slightly asymmetrically, to the 
Ca II K core. The data of April 9 have been dis- 
regarded. 



added in order to obtain a mean spectrum for each 
night. (4) Overall average spectra were computed over 
all nights. (5) Residuals given as {datai — mean) were 
computed. (6) Low-order curvature was removed from 
each residual spectra. This step was performed by ap- 
plying fits to the data sets with polynomials of order 2. 
In order to obtain the best trend line for the low-order 
curvature, data points in the proximity of the center 
line were excluded from the fit. (7) Next the resid- 
ual spectra is smoothed over 21 pixels. (8) Finally, 
the Mean Absolute Deviation is computed according 
to MAD = J2i=i 

ai — meanl with N denot- 
ing the number of spectra. As a further part of the 
process, all spectra were analyzed with respect to 7 A 
window. Trapezoidal integration was used to compute 
the integrated K residuals. Error bars were calculated 
from nightly variations. 

3 Results and Discussion 

The nightly fiux residuals in the K core are over-plotted 
in Fig. 3. The figure indicates fluctuations in the 
K core apparently due to star-planet i nteraction akin 
to the results previously obtained by IShkolnik et al.l 
(2003). The calculated Mean Absolute Deviation (MAD) 
is shown in Fig. 4. The average of all normalized 
nightly spectra (mean values) is depicted using dashed- 
lines indicating the MAD is confined, slightly asym- 
metrically, to the K core. A similar analysis has been 
pursued for the Ca II H line at 3968 A. The Ca II H 
line does not show any identifiable variation (see Figs. 
5 and 6). Previously, Ca II H fiuctuations in the line 



6 



Publications of the Astronomical Society of Australia 



-0.02 



-0.04 - 



3965 3966 3967 3968 3969 3970 3971 
Wovelength (A) 



Figure 5: Same as Fig. 3, but now for the Ca II H 
core of HD 179949. 



0.04 - 



0.02 - 



0.00 



-0.02 - 



-0.04 - 




3941 3942 3943 3944 3945 3946 3947 
Wovelength (A) 



< 0.01 - 



0.00 t 



0.04 - 



0.02 - 



0.00 



-0.02 - 



-0.04 - 




3965 3966 3967 3968 3969 3970 3971 
Wavelength (A) 



3941 3942 3943 3944 3945 3946 3947 
Wavelength (A) 



Figure 7: Top: Residual flux of the Al I A3944 A 

hne between April 5, 2006 and April 11, 2006 akin 
to Fig. 3 noting that the data of April 9 have 
again been disregarded. Bottom: Here we also dis- 
regarded the data of April 11, noting these data 
points either suffer from an unknown data artifact 
or are due to photospheric activity unrelated to 
star-planet interaction (see text). This signal has 
no counterpart in Ca II H and K. 



Figure 6: Same as Fig. 4, but now for the Ca II H 
core of HD 179949. 
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core were deduced to be about 2/3 of the streuRth o f 
the Ca II K core as stated by IShkolnik et al.l l|2003l '): 
see, e.g. iLinskv fc Avret3 (|l97C ) for background infor- 
mation on Ca II. We currently have no explanation for 
this discrepancy between the predicted and observed 
Ca II H emission. 

Another focus of our study is to investigate the be- 
havior of the Al I A3944 A hne (see Figs. 7 to 9), which 
is formed in the stellar photosphere. If the fluctuations 
in Ca II H, K are caused by star-planet interaction, as 
expected, rather than by an inherently stellar process, 
a significantly reduced level of variability in the stellar 
photosphere is expected to occur compared to the stel- 
lar chromosphere. The latter is a natural consequence 
of the strong density gradient between the stellar pho- 
tosphere and chrom osphere; see e.g. the V AL-C model 
for the Sun (G2 V) IVernazza et al.l (119811 ). which has 
a similar, albeit moderately lower surface temperature 
than HD 179949. Noting that Al I A3944 A is a photo- 
spheric line positioned between the Ca II H and K lines 
in the spectrum, it can thus be utilized for distinguish- 
ing between planet-related chromospheric activity and 
and planet-unrelated photospheric activity. 

The Al I line analysis is performed to isolate planet- 
induced chromospheric activity from possible photo- 
spheric fluctuations. Since the photospheric response 
is expected to be significantly less than the chromo- 
spheric signal, this is also a chec k on systematic errors 
in our data analysis. Previously. IShkolnik et al] (|2003l ) 
stated that no variation has been detected in 2 A win- 
dow of Al I analysis. We extend our analysis to a 7 A 
window to create results comparable to the Ca II H 
and K cores. It is found that the Al I line analysis in 
general does not indicate enhanced photospheric ac- 
tivity, even if enhanced chromospheric activity occurs. 
A notable exception, however, is the observation on 
April 11, 2006, where a strong Al I line variation is 
found without a chromospheric counterpart. Thus, if 
it is indicative of a real photospheric signal, it is ap- 
parently unassociated with star-planet interaction. 

Detected photospheric emission raises the level of 
the mean spectrum around the line core resulting in 
a greater residual flux (see Fig. 7). As noted in Sect. 
2, the residual flux is computed from (datUi — mean); 
therefore, the fluctuation increases because of an in- 
creased mean value. Integrated residuals between 3943.5 
A and 3944.5 A are plotted as a function of the plan- 
etary phase (see Fig. 8). The April 11 data, indicated 
by a square point, show no correlation with the ac- 
tivity enhancement detected in the K core of Ca II. 
If we ignore the April 11 data point regardless of the 
cause of the variability, and repeat our analysis, no 
signiflcant fluctuation and/or activity enhancement is 
revealed (see Fig. 7). The Mean Absolute Deviation 
(MAD) is subsequently also computed for the Al I line 
(Fig. 8) with no indication of activity enhancement 
found. 

To investigate the modulation of the PIE effect 
with the orbital period of the planet Porb, the Ca II K 
residuals are integrated and plotted as a function of 
orbital phase (see Fig. 10). The integrated residu- 
als are grouped by the orbital period of the planet. 
Each orbital coverage is distinguished by different sym- 
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Figure 8: Depiction of the integrated residual 
flux of the Al I A3944 A Hne of HD 179949. In 
the top figure, the data of April 9 have been dis- 
regarded, whereas in the bottom figure the data 
of both April 9 and April 11 have been omitted. 
Note that the data of April 11 are most likely due 
to photospheric activity unrelated to star-planet 
interaction. 
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Figure 9: Same as Fig. 4, but now for the Al I 
A3944 A line of HD 179949. Note that the data of 
both April 9 and April 11 have been disregarded 
(see text for details). 
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Table 4: Summary of Observations" 



Date 




UT 
(h) 


JD 




HJD 






m 


m 


<^(3) 


0(4) 


April 


6 2006 


11:31 


2453831 


.980 


2453831 


.980 


0, 


.2650 


0, 


.0137 


0, 


.1880 


0, 


.2750 


April 


6, 2006 


11:54 


2453831, 


.996 


2453831, 


.996 


0, 


.2702 


0, 


.0189 


0, 


.1931 


0, 


.2802 


April 


7, 2006 


11:17 


2453832, 


.970 


2453832 


.970 


0, 


.5852 


0, 


.3339 


0, 


.5082 


0, 


.5952 


April 


7, 2006 


11:40 


2453832, 


.986 


2453832 


.986 


0, 


.5904 


0, 


.3391 


0, 


.5133 


0, 


.6004 


April 


8, 2006 


11:30 


2453833 


.979 


2453833 


.979 


0, 


.9115 


0, 


.6602 


0, 


.8345 


0, 


.9216 


April 


8, 2006 


11:53 


2453833 


.995 


2453833 


.995 


0, 


.9167 


0, 


.6654 


0, 


.8396 


0, 


.9267 


April 


9, 2006 


11:13 


2453834, 


.967 


2453834 


.968 


0, 


.2311 


0, 


.9798 


0, 


.1540 


0, 


.2411 


April 


9, 2006 


11:36 


2453834, 


.983 


2453834, 


.983 


0, 


.2362 


0, 


.9849 


0, 


.1591 


0, 


.2462 


April 


9, 2006 


11:59 


2453834, 


.999 


2453835, 


.000 


0, 


.2414 


0, 


.9901 


0, 


.1643 


0, 


.2515 


April 


10, 2006 


11:59 


2453835, 


.999 


2453836 


.000 


0, 


.5647 


0, 


.3135 


0, 


.4877 





.5749 


April 


12, 2006 


11:31 


2453837, 


.980 


2453837, 


.980 


0, 


.2052 


0, 


.9540 


0, 


.1281 


0, 


.2154 


April 


12, 2006 


11:54 


2453837, 


.996 


2453837, 


.996 


0, 


.2103 


0, 


.9592 


0, 


.1333 


0, 


.2206 



"0(1), 0(2), 0(3), and 0(4) correspond to the observations obtained on 07/1998, 09/2000, 07/2002, and 
09/2003, respectively, (see Table 3) with 0(1) = 0ButiGr- 



b ols. The scale of the figure is adjusted compared 
to lShkolnik et all (|2003l ) by setting the minimal resid- 
ual flux to zero . The best-fit bright-spot models from 
IShkolnik et aP (|2003l 1 are also depicted; they imple- 
ment models of photospheric spatial inhomogeneities 
(spaced akin to Prot) in comparison to planetary effects 
(spaced akin to Porb)- The planetary positions for our 
observation run, calculated from the e phemerides of 
2002 provided bv IShkolnik et al.l (|2003l ). are depicted 
in Fig. 1. 

Figure 10 indicates enhanced activity during phases 
0.5 and 0.85 probably due to star-plan et interaction. 
This is consistent with previous results o f IShkolnik et al.l 
1^05, 2008) as t he data points are observed near bright- 
spot model fits. 'Shkolnik erill (|2005l . l2008f) suggested 
that the peak emission is found at a phase of about 
= 0.8 while stating that the fit to the 2001 and 
2002 data peaks at = 0.83 ± 0.04 with an amplitude 
of 0.027. The phase shifts in the HD 179949 system 
were first detected in previous work by E. Shkolnik and 
collaborators, when the observational data spanning 6 
years were analyzed. New fit functions suggested phase 
shifts of —0.07 between 2002 and 2005 data points, and 
-0.17 between 2003 and 2006 data points (Fig. 11). 

Currently, the cause(s) of these phase shifts are 
unknown, while noting that it may be possibly re- 
lated to the star-pla net interaction itse l f. As previ- 
ously pointed out (e.g. | McIvor et al.ll200d : IShkolnik et al.l 
120081 : ICohen et all [20091 '). the planet-induced activity 
leads the planet, and the amount of phase lead may 
permit to decipher the magnetic field geometry, in- 
cluding the existence of a Parker-type spiral. How- 
ever, considering that in HD 179949 the Parker spiral 
is relatively close to the star, a potentially more likely 
source of activity is magnetic reconnection to various 
flux-tubes at different longitudes and latitudes due to 
the tangled, vari able magnetic structures on the star 
(jCranmer &i Saaiii2008i V On the other hand, the phase 



shift may also be caused by a combination of differen- 
tial rotation and slow migrat ion of the ac tivity belt 
during a stellar dynamo cvcle [Lanzal (|2010l ). 

Note that we did not derive the ephemerides from 
our observa tional data, but used a v ailable 2002 ephemeri des 
provided bv lShkolnik et al.ll|2003l l. IShkolnik etaP (|2008l l 
reported phase shifts at an amount of —0.17 between 
the 2003 and 2006 observations; therefore, (almost) 
the same amount of shift may be applicable to our 
observations obtained in April 2006 at the McDonald 
Observatory. Figure 10 shows the results by including 
the —0.17 phase shift, which produces an appropri- 
ate correlation between our limited number of obser- 
vational data and the best-fit bright-spot models. The 
peak poi nt of the planet induced activity is estimated 
by Shko lnik et all (|2003l ^ to occur when the planet is 
about at the third quadrature (0 = 0.8). This primary 
estimation was based on best-fit bright-spot models 
over first available data points. However, the lack of 
a sufficient number of observational data and the de- 
tected phase shifts created a large uncertainty about 
the phase of the peak point. 

The latest four ephemerides of HD 179949 system 
are listed in Table 3. Calculated phases from the four 
different epochs are listed in Table 4. Table 5 shows 
the mean phase shifts (local time), noting that aver- 
ages were attained if there was more than one observa- 
tion per night. The inconsistency of calculated phases 
from different epochs may imply a change in mechan- 
ical properties in the system as function of time due 
to star- planet interaction. Figure 11 shows the change 
in phase as a function of time, where it shows a low- 
order polynomial trend. Data points are computed 
from A0 = — 0Buticr, where are the p hases given in 
Table 4 and 0Butior represent the data of lButler et al.l 
(|2006l '). The trend shown in Fig. 1 1 may be crucial for 
understanding the nature of the HD 179949 star-planet 
system; however, further discussions on this subject 
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0.04 f 



O First orbit 
A Second orbit 
□ Third orbit 





-0.5 0.0 0.5 1.0 

0.,. 



Figure 10: Top: Integrated residual flux of the 
Ca II K cores as a function of the orbital phase. 
Different symbols are used to distinguish the data 
points from the three consecutive orbits of HD 
179949b. All flux points are shifted vertically to 
be consistent with a minimum flux of zero. The 
ephemer ides for calcula. tinK the phases are those 
given bv lShkolnik et al.l (|2003h . The sohd line in- 
dic ates the best-fit bright -spot model as discussed 
by IShkolnik et al with a spot at a lati- 

tude of 30° and a stellar inclination of i = 87°. 
The dashed line represents the same model but for 
i = 83° . Note that the phase of the best-fit bright- 
spot models are shift ed by —0.17, an amount pre- 



viously suggested by IShkolnik et al.l ((2008) in the 
view of activity enhancement observations of 2003 
and 2006. Bottom: Same as top figure, but now 
for Ca II H. 




1998 1999 2000 2001 2002 2003 2004 
Yeors 



Figure 11: Phase shift difference relative to the 
observations obtained in July 1998. 



should be deferred until additional studies and obser- 
vations have become available. However, it is note- 
worthy that the existence of similar phase shifts have 
previously been identified in oth er star-planet system s, 
notably r Boo and HD 189733 l|Shkolnik et al.i l2008'l. 
possi bly indicativ e of stellar magnetic dynamo evolu- 
tion (|Lanzall2O10l '). 



Table 5: Mean Phase Shifts (Local Time)" 



Date 




0(2) 




m 


April 5, 2006 


0.2676 


0.0163 


0.1906 


0.2776 


April 6, 2006 


0.5878 


0.3365 


0.5108 


0.5978 


April 7, 2006 


0.9141 


0.6628 


0.8370 


0.9241 


April 8, 2006 


0.2362 


0.9850 


0.1591 


0.2463 


April 9, 2006 


0.5647 


0.3135 


0.4877 


0.5749 


April 11, 2006 


0.2078 


0.9566 


0.1307 


0.2180 



°0(1), (/)(2), (j){3), and 0(4) correspond to the 
observations obtained on 07/1998, 09/2000, 
07/2002, and 09/2003, respectively, (see Table 3) 

with 0(1) = 0Butlcr- 



4 Conclusions 

The aim of our study was to extend the observational 
database of planet induced stellar emission (PIE) due 
to the interaction of close-in exoplanets and their host 
stars by focusing on HD 179949. We obtained twelve 
observations during April 2006 at the 2.7 m telescope 
of the McDonald Observatory covering three planetary 
orbits. The PIE ph enomenon previously suggested by 
ICuntz et all (l2000l') w as first observed in HD 179949 
(jShkolnik et al.ll2003l ') concerning Ca II H and K. It is 
now kno wn to be present in five ot her star-planet sys- 
tems; see 'Cuntz fc ShkolnikI (2002) for an early com- 
parison between theory and observations. Besides Ca II 



10 



H and K, the PIE effect has also been identified in 

other sp ectral features in cluding the Ca II infr ared 

triplet fiSaar fc Cunt j2doih and in coronal X-rays (jKashvap 
l2008l : ISaar et al.ll2008l '). 

Our observations appear to confirm the activity en- 
han cement for the Ca II K core between phases 0.4 and 
0.9. IShkolnik et al.l (|2008h reported the disappearance 
of the PIE phenomenon during their 2003 and 2006 ob- 
servation runs. However, we have an insufficient num- 
ber of data to be able to confirm or refute this finding. 
Nevertheless, the two relevant data points obtained in 
our April 2006 observing run show the PIE effect to 
be present in the Ca II K core consistent with the 
previous results by Shkolnik and collaborators. Ad- 
ditionally, we obtained several data points showing 
the absence of star-planet interaction consistent with 
the previous predictions. This latter finding is consis- 
tent with the interpretation that the PIE effect is due 
to magnetic star-planet interaction presumably facili- 
tated along the Parker spiral of the stellar wind (e.g. 
ISaar et ahl 12004'). Recent magnetospheric simulations 
for the interaction between stellar and pl anetary mag- 
netic fields and outflows were given by ICohen et al.l 
l|201ll ) taking HD 189733 as an example, which point 
to the formation of associated Ca II line signatures, 
expected to occur in a stochastic manner. 

Regarding our observations, the Ca II H core anal- 
ysis did not reveal any significant fiuctuation, which 
is at odds with the prediction based on our Ca II K 
detection. Also, HD 179949 was observed at low alti- 
tude, hence at a relatively great air mass, during the 
observation run in April 2006, which is not suitable 
for detecting relativ ely weak act ivity enhancements. 
Previously, Shkolnik et al.l (120031) deduced the Ca II 
H core fluctuation to be about 2/3 of the intensity 
of the fluctuation observed in the K core. The Al I 
A3944 A line indicated no notable fluctuations during 
the time increments when the Ca II K emission en- 
hancements occurred. This is strong evidence that the 
observed emission is indeed due to the close-in giant 
planet, as expected. However, there is an exceptional 
Al I data point obtained on April 11 possibly revealing 
signiflcant emission in the stellar photosphere without 
a Ca II K chromospheric counterpart. Thus, the cause 
for this event is apparently inherent in the star, unre- 
lated to the existence of the planet. 
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Planet-Induced Emission Enhancements in HD 179949: 
Results from McDonald Observations 
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Abstract: We monitored the Ca II H and K lines of HD 179949 to observe and confirm previously 
identified planet induced emission (PIE) as an effect of star-planet interaction. We obtained high 
resolution spectra {R ~ 53,000) with a signal-to-noise ratio S/N > 50 in the Ca II H and K cores 
during 10 nights of observation at the McDonald Observatory. Wide band echelle spectra were taken 
using the 2.7 m telescope. Detailed statistical analysis of Ca II K revealed fiuctuations in the Ca II 
K core attr ibutable to pla n et ind uced chromospheric emission. This result is consistent with previous 
studies by IShkolnik et all l|2003l ). Additionally, we were able to confirm the reality and temporal 
evolution of the phase shift of the maximum of star-planet interaction previously found. However, no 
identifiable fiuctuations were detected in the Ca II H core. The Al I A3944 A line was also monitored to 
gauge if the expected activity enhancements are confined to the chromospheric layer. Our observations 
revealed some variability, which is apparently unassociated with planet induced activity. 



Keyvirords: planetary systems — radiation mechanisms: nonthermal 
mospheres — stars: individual (HD 179949) — stars: late-type 



stars: activity — stars: chro- 



1 Introduction 

Observational results indicate that the large major- 
ity of planets outside t he solar system is hosted by 
main-sequence stars fe.g. lButler et alll2006f ). which are 
known to possess an extended outer atmosphere con- 
sisting of a chromosphere, transition region and corona 
as well as a dynamic stellar wind region (e.g. iLinskvl 
Simon 1986) assoc iated with complex magnetic 
structure (e.g. Schrii ver fc ZwaanI [2000'). In various 
systems, consisting of a star and a close-in extrasolar 
giant planet (CEGP), interaction between the planet 
and stellar atmospheric structure has been identified 
that is often broadly classified as "planet-induced (stel- 
lar) emission" (PIE). The reali ty of this phenomenon 
has been proposed bv lCuntz et al. (2000) who also pro- 
vided a detailed ranking of the relative strength of this 
effect for the var ious star- planet sys tems known at the 
time. Moreover, I Cunt z et all (|200d ) explored the pos- 
sibility of tidal interaction and magnetic interaction 
with the latter distinguished from the former by two 
activity maxima per stellar circumference instead of 
one. 

Th e PIE effect was subsequently discovered bv lShkolnik 
(|2003l ') in regard to the HD 179949 system using high- 
resolution {R = 110, 000) and high signal-to-noise ratio 
spectra based on three observing runs at the Canada- 
France-Hawaii Telescope (CFHT). They also concluded 
that the PIE effect is expected to be magnetic in nature 
rather than caused by gravitational sta r-planet inter- 
action . Subsequent work was pursued bv lShkolnik et all 
l|2005l 'l, [Shkolnik et all (|2008l ') and others, resulting in 
the identification of the PIE effect in at least six sys- 
tems, which are: HD 209458, v And, r Boo, HD 179949, 
HD 189733, and HD 73526. 



The latter study also provides evidence of the on/off 
nature of the PIE effect, also previo usly found in the 
V And system l|Shkolnik et al.l [20081 ). which according 
to the authors is likely attributable to the changing 
stellar magnetic field structure thr oughout the stel- 
lar activity cycle (e.g. iLanzd [2OIOI ). An alternative 
or supplementary explanation of the significant time- 
dependent variations of the PIE phenome na refers to 
the fiare-type nature of th e interaction (e.g.lSa ar et al.l 
'20041: iMcIvor et al.! l2006^. Recently, a detailed time- 
dependent MHD simulation for the interaction of the 
stellar magnetic field and wind with the planetary mag- 
netosphere and outfiow was given bv lCohen et"all l|201ll ) 
for HD 189733 as an example. Their results show re- 
connection events occurring at specific planetary or- 
bital phases, causing mass loss from the planetary mag- 
netosphere able to generate hot spots on the stellar sur- 
face. The simulations also demonstrate that the sys- 
tem has sufficient energy for the hot spots to be visible 
in Ca II lines, which are expected to show stochastic 
behavior. 

There were also various attempts to identify the 
PIE effect in other spectral regi mes, including t he ra - 

tfiialtegime, infrared and X-rays. iBastian et all ([2OOOI ') 
conducted a search for radio emission for six extra- 
solar star-planet or brown dwarf systems that included 
four systems with the planet or BD closer than 0.1 AU 
from its host star. No detections were made, how- 
ever, which was attributed to physical or instrumen- 
tal reasons. Subsequent work mostly focused on the 
derivation of upper limits of eme rgent radio emissio n 
from known syste ms given bv. e.g. [ Farrell et al.l l|2003t ). 

IZarkal (|2007l ). and lGriefimeier et all l|2007f l. 

Anot her effort t o identify the PIE effect was under- 
taken bv lSaar fc Cuntz; (2001,) who executed a search 
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Tabic 1: Stellar and Planetary Parameters 



Parameter 


tiU 179949 


Keierencc 


Spectral Type 


F8.5 V 


a 


RA 


19 15 33.2278 


b 


DEC 


- 24 10 45.668 


b 


V 


6.25 ± ... 


b 


Tcff' (K) 


6202 ± 52 


Ribas at al. f2003) 


Af, (Mo) 


1.24 ± 0.10 


Tinnev et al. l'2001~) 


Distance (pc) 


27.0 ± 0.5 


b 


Age (Gyr) 


2.05 ± ... 


Donahue C 1993) 


R* (Rq) 


1.193 ± 0.030 


Ribas et al. f2003~) 


[Fe/H] 


0.226 ± 0.050 


Gonzalez & Laws ( 2007) 


ap (AU) 


0.0443 ± 0.0026 


Butler et al. (2006) 




0.022 ± 0.015 


Butler et al. (2006) 


Porb (days) 


3.092514 ± 0.000032 


Butler et al. (2006) 


MpS'mi (Mj) 


0.916 ± 0.076 


Butler et al. (2006) 



"data from SIMBAD; see |h ttp : / /simbad . u-strasbg . f rl 

adopted from the Hipparcos catalogue; see lESAI (|1997l) . 



for periodicities in the Ca II infrared triplet emission 
for a sample of stars with close-in giant planets; how- 
ever, no unequivocal identification was made. Addi- 
tionally, an effort to search for the PIE effect in the 
UV using the Goddard High Resolution Spectrograph 
(GHRS) on board of HST was rendered impossible as 
the GHRS was removed during HST Servicing Mission 
2. On the other hand, a statistical survey of the X-ray 
fluxes from stars with close-in planets has been fruitful 
as it found them enhanced by 30% to 40% on aver- 
age over typical fluxe s from similar stars w ith planets 
that are not close-in (jKashvap et al.ll2008l ). Direct X- 
ray o bservations of the HD 179949 system (jSaar et al.l 
120081 ) revealed a ~30% increase in the X-ray flux over 
quiescent levels coincident with the phase of the Ca II 
enhancements. Additionally, there is also a trend for 
the emission to be hotter at increased fluxes, confirmed 
by modeling, showing the enhancement at ~1 keV 
compared to ~0.4 keV for the background star. 

In our present study, we will again focus on the 
HD 179949 system. HD 179949 is an F8.5 V star 
at a distance of 27.0 ± 0.5 pc located in the south- 
ern hemisphere in the constellation Sagittarius with 
a magnitu de of V = 6.25. I ts mass is estimated as 
1.24 Mq (|Tinnev et all l200ll ) based on the interpo- 
lation of evolutionary t racks previously computed by 
iFuhrmann et all (|1997| . I1998D . The age of HD 179949 
has been estimated as 2.05 Gyr (see Table 1) based 
on the strength of the c hromospheric component of 
the Ca II emission, 

i?HK (|Noves et al.lll984l : iDonahud 
1 19931 ). Note that HD 179949 is a solar-type star of 
moderate age, implying that its inherent chromospheric 
emission and variability has considerably subsidecQ, 
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A 2nd orbit 
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Figure 1: Orbital position of the planet around 
its host star (HD 179949) during the observational 
run. Note that 4> — corresponds to the planetary 
position between the host star and Earth (sub- 
planetary point). The numbers indicate the dates 
in April 2006 (local time) when the observations 
were obtained. 



^See Guinan ct al. (2003) and iGiidell ||2007I 1. among 
other contributions to the "Sun in Time" project, for stud- 
ies on the decline of chromospheric emission on stellar evo- 
lutionary time scales. 



www.pu blish . csiro. au /journals/ pasa 



3 



0.8 -1 

- 

C 0.6 - 

CD 
> 

^04-1 
cr 

0.2 - 
0.0 L 



III 1^ 



3920 3930 3940 3950 3960 3970 39E 
Wavelength (A) 




3920 3930 3940 3950 3960 3970 398 
Wavelength(A) 



Figure 2: Single flat-fielded and blazed spectrum 
of HD 179949 taken on April 5, 2006 with an ex- 
posure time of 1200 s. We depict the Ca II regime 
(3915 ~ 3985 A) of the wide band echelle spectrum 
{top) and the associated S/N ratios {bottom). 



thus adding to the credibility of any PIE determina- 
tions. By employin g the Doppler velocity technique, 
iTinnev et alj (|200lD discovered as part of the Anglo- 
Australian Planet Search a "Hot Jupiter" hosted by 
HD 179949 located at a distance of about 0.044 AU 
(7.99 stellar radii). The planetary orbit has a very 
low eccentricity that is essentially consistent with the 
circular case (see Table 1), allowing us to ignore the or- 
bital ellipticity, if any, in the following. The planetary 
mass parameter is given as Mp sini = 0.916 ± 0.076 
Mj. 

Significant evidence that the activity enhancement 
in HD 179949 is planet-induced is based on the finding 
that the stellar rotation period Prot almost certainly 
disagees with the orbital period of the planet. There- 
fore, it would be difficult to explain emission on the 
star that essentially "follows" the CEGP without a 
planet-related cause. Unfortunately, HD 1799 49 docs 
not h ave a robustly measured Prot (see Wolf fc H armanec . 
120041 '). However, following iSaar et al] (|20oi )" and refer- 



ences therein, it can be argued based on the strength of 
the Ca II HK emission, the mean Mount Wilson index 
< S > as well as stellar evolution constraints that Prot 
lies between 7 and 10 d, thus being distinctly different 
from Porb =i 3.092 d (see Table 1). 

The PIE effect has various interesting astrophysi- 



cal applications as it, in principle, may allow the iden- 
tification and characterization of magnetic fields and 
dynamo activity in cxos olar Jupit er-type planets (e 
Id et al.l[20 04: Gricfimcie r et al.l20 04: Erkacv ct al. 200 
Zarkal2007i:)Lanza,200&lCohen et al..i200a : ,Schart,20l'^ 



Lanzall2011 : ICohen et al.ll2011^ . In addition, it may 
also be useful for pr obing the inner re gion of th e stel- 
lar winds (e.g. lPreus sc ct al. 2005; Mcl vor et al.i r2006). 
This is a further reason to continue monitoring the PIE 
effects in different spectral regimes and by using differ- 
ent type of instruments including the 2.7 m Harlan J. 
Smith Telescope at the McDonald Observatory. The 
latter allowed us to continue studying the Ca II H and 
K line cores in HD 179949 and to compare our findings 
with the previous results by Shkolink et al. (2003) and 
subsequent work. Our paper is structured as follows: 
In §2, we describe the employed methods including the 
data reduction. Our results and discussion are given 
in §3. In §4, we present our conclusions. 

2 Methods and Data Reduc- 
tion 

In order to study the PIE effect in the HD 179949 sys- 
tem, we obtained high resolution {R = A/ A A ~53,000) 
spectra using the McDonald Observatory's echelle spec- 
trograph mounted on the TK3 light path of the Coude 
focus of the 2.7 m (107 inch) Harlan J. Smith Tele- 
scope. We obtained observations for six (almost) con- 
secutive nights in April 2006, which allowed us to ac- 
quire data points spreading over nearly three planetary 
orbits (see Fig. 1 and Table 2 for details). Full optical 
bandwidth (3400 - 10,900 A) were recorded in at least 
67 orders by a Tektronix (2048 x 2048), CCD cam- 
era. The Ca II H and K lines were carefully centered 
in order 56. Furthermore, Th-Ar spectra were taken 
frequently (i.e., typically in less than 2 h intervals) to 
allow accurate wavelength calibrations. The signal-to- 
noise (S/N) ratio of the spectra exceeded 200 for the 
continuum, and it was close to 50 in the Ca II K core 
(3933 A), and about 80 in the Ca II H core (3968 A) 
(see Fig. 2). 

The data reduction was performed by using the 
IRAF standard packages IMRED, CCDRED and ECHELLE. 
The data reduction procedures can be summarized as 
follows: Combined biases were subtracted from all stel- 
lar, arc, and fiat-field exposures to remove the baseline 
noise. Scattered light was removed from all stellar and 
flat field exposures. All stellar and arc exposures were 
divided by combined flat fleld for normalization of the 
sensitivity levels of the CCD chip. Spectra were ex- 
tracted from the 2D images to 2D spectra. Heliocen- 
tric correction was applied to each spectrum to express 
the observation time considering that the observer's 
location does not coincide with the center of the Sun. 
Differential radial velocity corrections due to the rota- 
tion and orbital velocity of Earth were applied to each 
spectrum. The combined flat field was processed like 
a stellar image to obtain blaze functions for each ob- 
servation night. Blaze functions were used to remove 
from each flnal spectrum the low-order curvature as a 
result of the intensity variation along the orders. 
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Table 2: Details of Observations 



Number 


Date 


Time 


Exposure Time 


S/N 


Seeing 






(UT) 


(s) 




(arcs) 


1 


April 5, 2006 


11:31 


1200 


53 


2.6 


2 


April 5, 2006 


11:54 


1200 


56 


2.6 


3 


April 6, 2006 


11:17 


1200 


43 


4.1 


4 


April 6, 2006 


11:40 


1200 


47 


4.1 


5 


April 7, 2006 


11:30 


1200 


37 


5.0 


6 


April 7, 2006 


11:53 


1200 


45 


5.0 


7 


April 8, 2006 


11:13 


1200 


50 


2.3 


8 


April 8, 2006 


11:36 


1200 


51 


2.3 


9 


April 8, 2006 


11:59 


1200 


70 


2.3 


10 


April 9, 2006 


11:59 


1200 


44 


2.1 


11 


April 11, 2006 


11:31 


1200 


51 


1.7 


12 


April 11, 2006 


11:54 


1200 


66 


1.7 



Table 3: Ephemerides of HD 179949b 



Date 


JD 


Orbital Period 


Reference 






(d) 




07 / 1998 


2451001.510 


3.092514 


Butler et al. f2006)° 


09 / 2000 


2451793.980 


3.0925 


b 


07 / 2002 


2452479.823 


3.09285 


Shkolnik et al. C2003) 


09 / 2003 


2452894.110 


3.09246 


Shkolnik et al. C2003) 



"based on data obtained bv lTinnev et al. 1 (l200lh . 
'^^http : //www ■ exoplanet . eu 
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Figure 3: Residual flux of the Ca II K core of 
HD 179949 between April 5, 2006 and April 11, 
2006. The residual flux is computed from the av- 
erage mean and smoothed over 21 pixels. The re- 
sults are indicated by solid (April 5), dotted (April 
6), dashed (April 7), dash-dotted (April 8), and 
long-dashed (April 11) lines. Note that the data 
of April 9 have been disregarded. 
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Figure 4: The solid line indicates the Mean Ab- 
solute Deviation (MAD) of the Ca II K core of 
HD 179949 as observe d. The MAD va l ues ar e cal- 
culated as described bv lShkolnik et all (|2003[l . not- 
ing that MAD = X]i=i \datai — mean\ based 
on N spectra. The unit is the relative intensity 
attained as a fraction of the normalized flux. The 
dashed line represents the mean spectrum of the 
Ca II K core (scaled accordingly) to show that the 
activity is confined, slightly asymmetrically, to the 
Ca II K core. The data of April 9 have been dis- 
regarded. 



Figure 5: Same as Fig. 3, but now for the Ca II H 
core of HD 179949. 



Each reduced spectrum is carefully checked for po- 
tential contaminations. Trace amounts of sunlight con- 
tamination are detected in the last exposures of each 
observing night as HD 179949 is observed at low alti- 
tudes on the Eastern horizon just a few hours before 
dawn. These exposures were discarded to avoid any so- 
lar contamination of our spectra. Although the night 
of April 9 originally indicated a critical phase to de- 
tect star-planet interaction, this entire data set had to 
be discarded as the cloud coverage caused large noise 
level in the spectra. 



Ne xt we followed the procedure of IShkolnik et all 
l|2003l ) previously used to unveil the PIE effect for the 
HD 179949 system from optical spectra. The follow- 
ing steps were undertaken: (1) From each spectrum 
increments of 7 A centered at the Ca II H (3968 A) 
and K (3933 A) cores, respectively, were extracted. 
(2) To normalize each spectrum, the end points of 
each extracted spect ra were set to 1. As stated by 
IShkolnik etaP (|2003l i. this is the best way to normal- 
ize the spectra in the wavelength regimes, such as H 
and K lines, where there is no clear continuum. (3) 
The normalized spectra were grouped by date and co- 
added in order to obtain a mean spectrum for each 
night. (4) Overall average spectra were computed over 
all nights. (5) Residuals given as {datUi — mean) were 
computed. (6) Low-order curvature was removed from 
each residual spectra. This step was performed by ap- 
plying fits to the data sets with polynomimals of or- 
der 2. In order to obtain the best trend line for the 
low-order curvature, data points in the proximity of 
the center line were excluded from the fit. (7) Next 
the residual spectra is smoothed over 21 pixels. (8) 
Finally, the Mean Absolute Deviation is computed ac- 
cording to MAD = iV"^ J^'Li \dato,r 

— mean] with A'^ 

denoting the number of spectra. As a further part 
of the process, all spectra were analyzed with respect 
to 7 A window. Trapezoidal integration was used to 
compute the integrated K residuals. Error bars were 
calculated from nightly variations. 
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Figure 6: Same as Fig. 4, but now for the Ca II H 
core of HD 179949. 
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Figure 7: Top; Residual flux of the Al I A3944 A 
line between April 5, 2006 and April 11, 2006 akin 
to Fig. 3 noting that the data of April 9 have 
again been disregarded. Bottom: Here we also dis- 
regarded the data of April 11, noting these data 
points either suffer from an unknown data artifact 
or are due to photospheric activity unrelated to 
star-planet interaction (see text). This signal has 
no counterpart in Ca II H and K. 



3 Results and Discussion 

The nightly flux residuals in the K core are over-plotted 
in Fig. 3. The figure indicates fluctuations in the 
K core apparently due to star-planet i nteraction akin 
to the results previously obtained by IShkolnik et al.l 
(|2003l '). The calculated Mean Absolute Deviation (MAD) 
is shown in Fig. 4. The average of all normalized 
nightly spectra (mean values) is depicted using dashed- 
lines indicating the MAD is confined, slightly asym- 
metrically, to the K core. A similar analysis has been 
pursued for the Ca II H hue at 3968 A. The Ca II H line 
does not show any identifiable variation. Previously, 
Ca II H fluctuations in the line core were deduced to be 
ab out 2/3 of the st r ength of the Ca I I K core as stated 
bv IShkolnik et all (|2003l '): see, e.g. iLinskv fc AvrettI 
(|197(]I ') for background information on Ca II. We cur- 
rently have no explanation for this discrepancy be- 
tween the predicted and observed Ca II H emission. 

Another focus of our study is to investigate the be- 
havior of the Al I A3944 A line (see Figs. 7 to 9), which 
is formed in the stellar photosphere. If the fluctuations 
in Ca II H, K are caused by star-planet interaction, as 
expected, rather than by an inherently stellar process, 
a significantly reduced level of variability in the stellar 
photosphere is expected to occur compared to the stel- 
lar chromosphere. The latter is a natural consequence 
of the strong density gradient between the stellar pho- 
tosphere and chrom osphere; see e.g. the V AL-C model 
for the Sun (G2 V) IVernazza et al.l (ligsil ). which has 
a similar, albeit moderately lower surface temperature 
than HD 179949. Noting that Al I A3944 A is a photo- 
spheric line positioned between the Ca II H and K lines 
in the spectrum, it can thus be utilized for distinguish- 
ing between planet-related chromospheric activity and 
and planet-unrelated photospheric activity. 

The Al I line analysis is performed to isolate planet- 
induced chromospheric activity from possible photo- 
spheric fluctuations. Since the photospheric response 
is expected to be signiflcantly less than the chromo- 
spheric signal, this is also a chec k on systematic errors 
in our data analysis. Previously IShkolnik et all (|2003l ) 
stated that no variation has been detected in 2 A win- 
dow of Al I analysis. We extend our analysis to a 7 A 
window to create results comparable to the Ca II H 
and K cores. It is found that the Al I line analysis in 
general does not indicate enhanced photospheric ac- 
tivity, even if enhanced chromospheric activity occurs. 
A notable exception, however, is the observation on 
April 11, 2006, where a strong Al I line variation is 
found without a chromospheric counterpart. Thus, if 
it is indicative of a real photospheric signal, it is ap- 
parently unassociated with star-planet interaction. 

Detected photospheric emission raises the level of 
the mean spectrum around the line core resulting in 
a greater residual flux (see Fig. 7). As noted in Sect. 
2, the residual flux is computed from {datUi — mean); 
therefore, the fluctuation increases because of an in- 
creased mean value. Integrated residuals between 3943.5 
A and 3944.5 A are plotted as a function of the plan- 
etary phase (see Fig. 8). The April 11 data, indicated 
by a square point, show no correlation with the ac- 
tivity enhancement detected in the K core of Ca II. 
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Figure 8: Depiction of the integrated residual 
flux of the Al I A3944 A Hue of HD 179949. In 
the top figure, the data of April 9 have been dis- 
regarded, whereas in the bottom figure the data 
of both April 9 and April 11 have been omitted. 
Note that the data of April 11 are most likely due 
to photospheric activity unrelated to star-planet 
interaction. 



If we ignore the April 11 data point regardless of the 
cause of the variability, and repeat our analysis, no 
significant fluctuation and/or activity enhancement is 
revealed (see Fig. 7). The Mean Absolute Deviation 
(MAD) is subsequently also computed for the Al I line 
(Fig. 8) with no indication of activity enhancement 
found. 

To investigate the modulation of the PIE effect 
with the orbital period of the planet Porb, the Ca II K 
residuals are integrated and plotted as a function of 
orbital phase (see Fig. 10). The integrated residu- 
als are grouped by the orbital period of the planet. 
Each orbital coverage is distinguished by different sym- 
bo ls. The scale of the figure is adjusted compared 
to lShkolnik et all l|2003l ) by setting the minimal resid- 
ual flux to zero . The best-fit bright-spot models from 
IShkolnik et al.l (|2003l l are also depicted; they imple- 
ment models of photospheric spatial inhomogeneities 
(spaced akin to Prot) in comparison to planetary effects 
(spaced akin to Porb). The planetary positions for our 
observation run, calculated from the e phemerides of 
2002 provided by IShkolnik et al.l (|2003l l. are depicted 
in Fig. 1. 



Figure 9: Same as Fig. 4, but now for the Al I 
A3944 A line of HD 179949. Note that the data of 
both April 9 and April 11 have been disregarded 
(see text for details). 



Figure 10 indicates enhanced activity during phases 
0.5 and 0.85 probably due to star-plan et interaction. 
This is cons istent with previous results o f IShkolnik et al.l 
(|2005l . l2008l ) as t he data points are observed near bright- 
spot model fits. IShkolnik eral . (2005, 2008) suggested 
that the peak emission is found at a phase of about 
(j> = 0.8 while stating that the fit to the 2001 and 
2002 data peaks a,t (j) = 0.83 ± 0.04 with an amplitude 
of 0.027. The phase shifts in the HD 179949 system 
were first detected in previous work by E. Shkolnik and 
collaborators, when the observational data spanning 6 
years were analyzed. New fit functions suggested phase 
shifts of —0.07 between 2002 and 2005 data points, and 
-0.17 between 2003 and 2006 data points (Fig. 11). 

Currently, the cause(s) of these phase shifts are 
unknown, while noting that it may be possibly re- 
lated to the star-pla net interacti on itself. As prevl 



ously pointed out (e.g. | McIvor etl il. 2006; S hkolnik et al.1 
l2008l : ICohen et al.ll2009f ). the planet-induced activity 
leads the planet, and the amount of phase lead may 
permit to decipher the magnetic field geometry, in- 
cluding the existence of a Parker-type spiral. How- 
ever, considering that in HD 179949 the Parker spiral 
is relatively close to the star, a potentially more likely 
source of activity is magnetic reconnection to various 
flux-tubes at different longitudes and latitudes due to 
the tangled, variable m agnetic structures on the star 
(jCranmer fc Saarll2008f ). On the other hand, the phase 
shift may also be caused by a combination of differen- 
tial rotation and slow migrat ion of the ac tivity belt 
during a stellar dynamo cvcle [banzai (|2010l b 

Note that we did not derive the ephemerides from 
our observational data, but used available 2002 ephemerides 
provided by Shkolnik et al. (2003). Shkolnik ct al. (2008) 
reported phase shifts at an amount of —0.17 between 
the 2003 and 2006 observations; therefore, (almost) 
the same amount of shift may be applicable to our 
observations obtained in April 2006 at the McDonald 
Observatory. Figure 10 shows the results by including 
the —0.17 phase shift, which produces an appropri- 
ate correlation between our limited number of obser- 
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Table 4: Summary of Observations" 



Date 




UT 
(h) 


JD 




HJD 






m 




<^(3) 


0(4) 


April 


6 2006 


11:31 


2453831 


980 


2453831 


980 





2650 


0.0137 





1880 


0.2750 


April 


6, 2006 


11:54 


2453831 


996 


2453831 


996 





2702 


0.0189 





1931 


0.2802 


April 


7, 2006 


11:17 


2453832 


970 


2453832 


970 





5852 


0.3339 





5082 


0.5952 


April 


7, 2006 


11:40 


2453832 


986 


2453832 


986 





5904 


0.3391 





5133 


0.6004 


April 


8, 2006 


11:30 


2453833 


979 


2453833 


979 





9115 


0.6602 





8345 


0.9216 


April 


8, 2006 


11:53 


2453833 


995 


2453833 


995 





9167 


0.6654 





8396 


0.9267 


April 


9, 2006 


11:13 


2453834 


967 


2453834 


968 





2311 


0.9798 





1540 


0.2411 


April 


9, 2006 


11:36 


2453834 


983 


2453834 


983 





2362 


0.9849 





1591 


0.2462 


April 


9, 2006 


11:59 


2453834 


999 


2453835 


000 





2414 


0.9901 





1643 


0.2515 


April 


10, 2006 


11:59 


2453835 


999 


2453836 


000 





5647 


0.3135 





4877 


0.5749 


April 


12, 2006 


11:31 


2453837 


980 


2453837 


980 





2052 


0.9540 





1281 


0.2154 


April 


12, 2006 


11:54 


2453837 


996 


2453837 


996 





2103 


0.9592 





1333 


0.2206 



"0(1), 0(2), 0(3), and 0(4) correspond to the observations obtained on 07/1998, 09/2000, 07/2002, and 
09/2003, respectively, (see Table 3) with 0(1) = 0ButiGr- 



vational data and the best-fit bright-spot models. The 
pe ak point of the planet induced activity is estimated 
bv lShkolnik et all (|2003t) to occur when the planet is 
about at the third quadrature (0 = 0.8). This primary 
estimation was based on best-fit bright-spot models 
over first available data points. However, the lack of 
a sufficient number of observational data and the de- 
tected phase shifts created a large uncertainty about 
the phase of the peak point. 

The latest four ephemerides of HD 179949 system 
are listed in Table 3. Calculated phases from the four 
different epochs are listed in Table 4. Table 5 shows 
the mean phase shifts (local time), noting that aver- 
ages were attained if there was more than one observa- 
tion per night. The inconsistency of calculated phases 
from different epochs may imply a change in mechan- 
ical properties in the system as function of time due 
to star- planet interaction. Figure 11 shows the change 
in phase as a function of time, where it shows a low- 
order polynomial trend. Data points are computed 
from A0 = — 0ButiGr, where are the p hases given in 
Table 4 and 0Butior represent the data of lButler et al.l 
(|2006l ). The trend shown in Fig. 11 may be crucial for 
understanding the nature of the HD 179949 star-planet 
system; however, further discussions on this subject 
should be deferred until additional studies and obser- 
vations have become available. However, it is note- 
worthy that the existence of similar phase shifts have 
previously been identified in oth er star-planet system s, 
notably r Boo and HD 189733 (|Shkolnik et al.i 120081 ') . 
possibly indicative of stellar magnetic dynamo evolu- 
tion (,Lanza,201Q '). 

4 Conclusions 

The aim of our study was to extend the observational 
database of planet induced stellar emission (PIE) due 
to the interaction of close-in exoplanets and their host 



Table 5: Mean Phase Shifts (Local Time)" 



Date 




0(2) 


0(3) 


0(4) 


April 5, 2006 


0.2676 


0.0163 


0.1906 


0.2776 


April 6, 2006 


0.5878 


0.3365 


0.5108 


0.5978 


April 7, 2006 


0.9141 


0.6628 


0.8370 


0.9241 


April 8, 2006 


0.2362 


0.9850 


0.1591 


0.2463 


April 9, 2006 


0.5647 


0.3135 


0.4877 


0.5749 


April 11, 2006 


0.2078 


0.9566 


0.1307 


0.2180 



"0(1), 0(2), 0(3), and 0(4) correspond to the 
observations obtained on 07/1998, 09/2000, 
07/2002, and 09/2003, respectively, (see Table 3) 

with 0(1) EE 0Butlcr- 



stars by focusing on HD 179949. We obtained twelve 
observations during April 2006 at the 2.7 m telescope 
of the McDonald Observatory covering three planetary 
orbits. The PIE ph enomenon previously suggested by 
LCuntz et al.' (2000) was first observed in HD 179949 
(Shkolnik ot al. 2003i) concerning Ca II H and K. It is 
now kno wn to be present in f ive ot her star-planet sys- 
tems; see ICuntz fc ShkolnikI (|2002l ) for an early com- 
parison between theory and observations. Besides Ca II 
H and K, the PIE effect has also been identified in 

other s pectral features inc luding the Ca II infr ared 

triplet (Saar & Cuntz 200i|) and in coronal X-rays (|Kashvap et al.l 
r200& : .Saar et al, 2008). 

Our observations appear to confirm the activity en- 
han cement for the Ca II K core between phases 0.4 and 
0.9. IShkolnik et al.l (|2008l ) reported the disappearance 
of the PIE phenomenon during their 2003 and 2006 ob- 
servation runs. However, we have an insufficient num- 
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Figure 10: Top: Integrated residual flux of the 
Ca II K cores as a function of the orbital phase. 
Different symbols are used to distinguish the data 
points from the three consecutive orbits of HD 
179949b. All flux points are shifted vertically to 
be consistent with a minimum flux of zero. The 
ephemer ides for calcula. tinK the phases are those 
given bv lShkolnik et al.l (|2003h . The sohd line in- 
dic ates the best-fit bright -spot model as discussed 
by IShkolnik et al with a spot at a lati- 

tude of 30° and a stellar inclination of i = 87°. 
The dashed line represents the same model but for 
i = 83° . Note that the phase of the best-fit bright- 
spot models are shift ed by —0.17, an amount pre- 



viously suggested by IShkolnik et al.l ( j2008) in the 
view of activity enhancement observations of 2003 
and 2006. Bottom: Same as top figure, but now 
for Ca II H. 
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Figure 11: Phase shift difference relative to the 
observations obtained in July 1998. 



ber of data to be able to confirm or refute this finding. 
Nevertheless, the two relevant data points obtained in 
our April 2006 observing run show the PIE efi'ect to 
be present in the Ca II K core consistent with the 
previous results by Shkolnik and collaborators. Ad- 
ditionally, we obtained several data points showing 
the absence of star-planet interaction consistent with 
the previous predictions. This latter finding is consis- 
tent with the interpretation that the PIE effect is due 
to magnetic star-planet interaction presumably facili- 
tated along the Parker spiral of the stellar wind (e.g. 
ISaar et aP 2004). Recent magnetospheric simulations 
for the interation between stellar and pl anetary mag- 
netic fields and outflows were given by ICohen et al.l 
(|201lh taking HD 189733 as an example, which point 
to the formation of associated Ca II line signatures, 
expected to occur in a stochastic manner. 

Regarding our observations, the Ca II H core anal- 
ysis did not reveal any significant fluctuation, which 
is at odds with the prediction based on our Ca II K 
detection. Also, HD 179949 was observed at low alti- 
tude, hence at a relatively great air mass, during the 
observation run in April 2006, which is not suitable 
for detecting rela tively we ak activity enhancements. 
Previously, Shkol nik et all ((200a 1 deduced the Ca II 
H core fluctuation to be about 2/3 of the intensity 
of the fluctuation observed in the K core. The Al I 
A3944 A line indicated no notable fluctuations during 
the time increments when the Ca II K emission en- 
hancements occurred. This is strong evidence that the 
observed emission is indeed due to the close-in giant 
planet, as expected. However, there is an exceptional 
Al I data point obtained on April 11 possibly revealing 
signiflcant emission in the stellar photosphere without 
a Ca II K chromospheric counterpart. Thus, the cause 
for this event is apparently inherent in the star, unre- 
lated to the existence of the planet. 
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